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Abstract 

We investigate low temperature properties of a random Ising model with + J 

and —aJ{a ^ 1) bonds in two dimensions using a cluster heat bath method. 

It is found that the Binder parameters gi for different sizes of the lattice 

come together at almost the same temperature implying the occurrence of the 

spin glass(SG) phase transition. From results of finite size scaling analyses, 

we suggest that the SG phase really occurs at low temperatures which is 

characterized by a power law decay of spin correlations. 
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Spin glasses have attracted great challenge for computational physics in these two 
decades. It is widely believed now in the bond-random Ising model that spin glass (SG) 
transitions occur at a finite, non-zero temperature Tc 7^ in three dimensions(3D) and 
at zero temperature Tc = in two dimensions(2D) @-§|. Recently, the present authors 
reexamined the SG phase transition of the ±J Ising model on a square lattice of L x L by 
means of an exchange Monte Carlo method and found that the Binder parameters ql 
for L < 16 intersect at T 7^ 0. They also found that better finite-size scaling(FSS) fits of the 
spin glass susceptibility xsg are obtained when Tc 7^ 0. These results imply the occurrence 
of the SG phase transition at Tc 7^ 0. If so, it is quite interesting, because it disproves 
the belief of Tc = 0. However, there remain two problems which should be considered to 
suggest Tc 7^ in 2D. One is that ql for a smaller lattice almost saturates below a rather 
high temperature [0 and its saturation value slightly increases with L 0. Therefore, it is 
difficult to see whether the intersection of ql for smaller L suggests the presence of the SG 
phase at Tc 7^ or is merely due to a finite size effect. The other is that it is still open 
whether or not the model really exhibits the nature of the SG phase at T < Tc, because the 
estimated transition temperature Tc is a slightly lower than the lowest temperature which 
is reached in the simulation. The problems would be solved, if we study the same model 
on bigger lattices at lower temperatures. The saturation oi ql at rather high temperatures, 
however, may be removed, if we treat an asymmetric random Ising model with +J and 
—aJ{a 7^ 1) bonds, because the energy gap of 2|1 — a| J in that model between the ground 



state and the lowest excitation state is much smaller than that of 4J in the ±J model 11 



and, if the lattice is rather small, we may study equilibrium properties at any temperature 



using a cluster heat bath(CHB) method 1|T2|,|T3 



In this Letter, we investigate low temperature properties of the asymmetric random Ising 
model on the square lattice of T x L{L < 18) using the CHB method. In fact, does not 
saturate down to a very low temperature. We find that, as the temperature decreases, 
qls for different L meet at almost the same temperature and then increase together. This 
property rather resembles that of the ±J model in 3D in which the SG phase transition 



occurs at Tc 7^ 0. We make the FSS analyses and find that ql and xsg for different L scale 
well using a finite, non-zero value of and that the distribution functions Pl{Q) of the 
spin overlap Q scale at all temperatures below Tc. Thus we suggest that in this model the 
SG phase occurs at low temperatures which is characterized by a power law decay of the 
spin correlations. The properties for T > found here are very similar to those of the ± J 
model in 2D [^. We believe, hence, that the SG phase transition also occurs at Tc 7^ in 
the ±J model in 2D. 

We start with an Ising model on a square lattice L x L described by the Hamiltonian 

H = - -Jij^^i^^ji (1) 

<i,j> 

where crj(= ±1) are Ising spins and < ij > runs all nearest neighbor pairs. Distributions of 
bonds for different samples are given as 

PiJr,) = lmj-J) + 6iJ,,+aJ)]. (2) 

If the lattice is rather small and a free boundary condition is used at least for one direction, 
we may readily obtain equilibrium spin configurations at any temperature T by using the 
CHB method, because the cluster defined in Ref. 13 can be chosen as the lattice itself. That 
is, the cluster is composed of L layers with L spins and exchange fields from the outside are 
absent. We briefly note the method For every sample, the weight functions Fi{{af^}) 
can be uniquely determined by using Eq.(7) in Ref. 13, because hf^ = 0. Once the set of 
these functions {T)({crf "*})} is obtained, the spin configurations of individual layers can be 
determined successively from the layer (L) to the layer (1) by using Eq.(8) or Eq.(ll) in Ref. 
13 and random numbers. Thus one of the spin configurations of the lattice is generated. 
Repeating this procedure, we may generate any number of the spin configurations with 
the aid of {Fi{{af^})}. These spin configurations are independent with each other and in 



accordance with the Boltzmann's weight |T3[. For each of the samples, about M = 200 spin 
configurations are generated ||T5|. We calculate, as well as usual magnetic quantities, an 
overlap function of the spins between the different spin configurations, Pi^\Q), for every 
sample: 



M M 



(Q) = Mij^) T.T.KQ-Q':i) (3) 
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with = (1/^) Z^ili where o"!"'^'* is the ith spin of the nth spin con- 

figuration for the fcth sample. The overlap function Pl{Q) of the system is given as 
Pl{Q) = / Ns) J2k=i Pl'\Q) ^ where Ns is the number of the samples. Once Pl{Q) is 
determined, we may obtain various SG quantities. The nth moment of the spin overlap is 
defined as: 

[<g">] = J^' Q'^PL{Q)dQ, (4) 

where < • • • > and [■ ■ ■] mean the spin configuration(thermal) average and the sample 
average, respectively. The SG susceptibility xsg is determined from 

XSG = N[< >], (5) 

and the Binder parameter ql from 

g^ = i3-[<Q^>]/[<Q'>f)/2. (6) 

We have performed this CHB simulation of the model of Eq.(l) with a = 0.8 for L < 18. 
The numbers of the samples are Ng = 4000 for L < 14 and Ng = 1000 for L = 16 and 18. 

Figure 1 shows Pl{Q) for different L. For every size L, the shape is symmetric with 
respect to Q = and the peaks at Q ~ ±1 become steeper as the temperature decreases. 
Figure 2 shows plots of gi against T. As the temperature decreases, gis for different L meet 
at almost the same temperature of T ~ 0.2J and then increase together. This behavior is 
quite similar to that of g^ of the ± J model in 3D in which the SG phase transition occurs 
at T, ^ i. 

We examine the results in more detail using the FSS analyses. First, we perform the 
FSS plots of gL and xsg to estimate the value of Tc. If the SG transition occurs at Tc, gi 
and Xsg scale as 

g^ = GieL'/n, (7) 
XsG = L'-^X{eL'/n, (8) 



where e = (T — Tc) /J, z/ is the exponent of the correlation length, t] is the exponent which 
describes the decay of the spin correlation at T = T^, and G and X are some scaling 
functions. Having assumed Tc ^ 0.1 9 J, we could obtain good scaling plots for T > Tc. 
Typical examples are shown in Figs. 3(a) and 3(b). Of course, the values of u estimated 
from both the scaling plots are almost the same. We also examined the other possibility of 
Tc = 0. As for XSG, we could scale the data only in the neighborhood of T = using 77 ~ 
and u ~ 2.6. As for ql, on the other hand, we could never scale the data using any plausible 
value of u. The scaling plots for Tc = are shown in Figs. 4(a) and 4(b). These results 
clearly reveal that, if a conventional phase transition occurs, the transition temperature is 
Tc ~ 0.19 J, not Tc ~ 0. However, the data for T < Tc deviate from the scaling plots. 

Next we examine Pl{Q) itself to see whether or not the SG phase is realized below 
T < Tc- If the phase transition occurs at T = T^, Pl{Q)^s for different L will scale as 

PL{q) = L'^/'P{qL^/') at T = T,. (9) 

Since Pl{Q) for smaller L has a considerable weight at Q = ±1 and, for L < 16, the peak 
height rather decreases with increasing L, we could not scale all the data over the entire 
range of Q. These difficulties may, however, come from finite size effects. In fact, the weights 
aX Q = ±1 become smaller as L increases and, for T < 0.2J, the peak height for L = 18 
becomes larger than that for L = 16. If we overlook the discrepancy around the peak of 
Pl{Q), the data scale for L > 10 at T ~ Tc by using the same value of 77 ~ 0.14 for Xsg, 
which is shown in Fig. 5(a). However, the data for T < Tc can also be scaled by using 
a smaller value of rj as shown in Fig. 5(b) |]16[. A natural interpretation of this result is 
that the model is close to criticality at all temperatures below Tc like the XY ferromagnet 



in 2D |[T^. This picture is, of course, compatible with the fact that and xsg scale only 
for T >Tc. We suggest, hence, that the SG phase really occurs below Tc ~ 0.19 J which is 
characterized by a power law decay of the spin correlations. 

Our present result is in agreement with our previous result of Tc 7^ in the ±J model 0. 
Especially, the value of z/ ~ 0.18 in the ±J model is in good agreement with that obtained 



in the present model [|T8|- Thus we predict that the occurrence of the SG phase at Tc 7^ 



is the common nature of 2D random Ising models with a discrete distribution of bonds [jT9[ . 
Our prediction of Tc 7^ appears incompatible with the previous belief of Tc = 0. We 
think, however, that these are not necessarily incompatible, because the previous authors 
have only concluded that their data are not incompatible with the prediction that = 0, 
and have not ruled out the possibility of such low Tc as estimated here. The thing that was 
certainly suggested by the previous studies is that, at T = 0, no long-range order exists and 
the spin correlation decays according to the power law This is compatible with our 

result with Tc ^ 0. Of course, further studies are necessary to confirm the prediction of 
Tc 7^ 0. We believe that the present result will stimulate not only the computational physics 
but also experimental studies, because the bond distributions of real SG materials will be 
asymmetric. 

One of the authors (TS) wishes to thank Professor H. Takayama, Professor H. Kawamura, 
Dr. Y. Ozeki, Dr. K. Hukushima and Dr. H. Yoshino for valuable discussions. The 
simulations were made partly on FACOM VPP500 at the Institute for Solid State Physics 
in University of Tokyo. 
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FIGURES 

FIG. 1. Pl(Q) versus Q at (a) T = 0.3J, (b) T = 0.2J and (c) T = O.IJ in the random Ising 
model with a = 0.8 on the square lattice oi L x L. 

FIG. 2. Temperature dependences of gL of the random Ising model with a = 0.8. 

FIG. 3. Scaling plots of (a) gi and (b) XSG, assuming Tc = 0.19 J and e = (T — Tc)/J. 

FIG. 4. Scaling plots of (a) gi and (b) xsG-, assuming = and e = (T — T^)/ J ■ 

FIG. 5. Scaling plots of Pl{Q) at (a) T = 0.2J and (b) T = O.IJ. 
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